Robust approach to f(R) gravity 
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We consider metric f(R) theories of gravity without mapping them to their scalar-tensor coun- 
terpart, but using the Ricci scalar itself as an "extra" degree of freedom. This approach avoids 
then the introduction of a scalar- field potential that might be ill defined (not single valued). In 
order to explicitly show the usefulness of this method, we focus on static and spherically symmetric 
spacetimes and deal with the recent controversy about the existence of extended relativistic objects 
in certain class of f(R) models. 
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I. INTRODUCTION 

Modified f(R) theories of gravity are metric theories 
which postulate an a priori arbitrary function of the Ricci 
scalar R as the Lagrangian density. These theories have 
been proposed to explain the late time accelerated expan- 
sion of the Universe as well as a mechanism to produce 
inflation in the early Universe in terms of geometry in- 
stead of introducing any dark energy ideas or a scalar in- 
flaton Thanks to the above properties, these theories 
have become one of the most popular alternative theories 
of gravity. In recent years, several specific f(R) models 
have been analyzed in different settings (see Ref. 0, 
for a review). Although the early models were ruled 
out for failing several tests (like the Solar System ex- 
periments) new proposals were put forward to overcome 
such drawbacks. Nevertheless, a considerable amount of 
analysis and observational confrontation is still required 
in order to compare the preliminary successes of such 
theories with the great achievements of general relativ- 
ity (GR). In particular, the models that were claimed 
to pass several cosmological and local tests have been 
scrutinized in the strong gravity regime only recently. In 
fact, the first studies concerning their ability to describe 
relativistic extended objects, like neutron stars, seem to 
give contradictory results. Using a f(R) model proposed 
by Starobinsky [3] , Kobayashi and Maeda Q found that 
such relativistic objects were difficult to construct since a 
curvature singularity developed within the object. Later, 
this issue was reanalyzed by several authors [H, , who 
found that singularity-free relativistic objects can indeed 
be constructed, arguing that the conclusion reached in 
[H was due to a bad choice of the matter model (i. e. 
the equation of state, hereafter EOS) Jj|, while in @ it 
was claimed that a "chameleon" is the responsible for the 
existence of such "stars," regardless of the EOS. 

A common feature to all of the aforementioned 
works 04^] is that the analysis was performed by map- 
ping the Starobinsky model to a scalar-tensor theory of 



gravity (STT) (c/. [7]). Under this mapping, the Ricci 
scalar has a behavior of the sort R ~ l/(x — Xo)> where 
X '■= diif is the scalar-field associated with the corre- 
sponding STT and xa — const. The key point is to de- 
termine if the dynamics of x leads it or not to the value 
X = Xo within the spacetime generated by the relativistic 
object. Irrespective of the different results and confusing 
explanations obtained in [U-[6j], we will argue that their 
conclusions are rather questionable due to the fact that 
the mapping to the STT is ill defined. To be more spe- 
cific, the scalar-field potential used to study the dynamics 
of x is n °t single valued and possesses pathological fea- 
tures. Since similar kind of singularities were also found 
in the cosmological setting Q , it is then worrisome that 
the ill-defined potential play such a crucial role in those 
analyses (c/. Refs. for further criticisms. See also 

and references therein for a more detailed discussion 
on cosmological singularities in f(R) theories). 

The goal of this communication is threefold: 1) to 
propose a straightforward and robust approach which 
consists in recasting the field equations in a more suit- 
able way without mapping the original f(R) theory to 
any scalar-tensor counterpart. This method dispenses us 
from dealing with ill-defined quantities that might arise 
when performing such transformation; 2) to reanalyze 
carefully the issue about the existence of relativistic ex- 
tended objects using our approach; 3) to stress that the 
analysis of f(R) theories based on the STT analogue with 
ill-defined quantities is not trustworthy, and that in cases 
where the STT approach is well defined (this depends on 
the specific f(R) model) it turns out to be rather convo- 
luted and not very insightful. Since f(R) theories are still 
under close examination, a sound approach is required to 
treat them appropriately. This is the first step in that 
direction. 



II. f(R) THEORIES OF GRAVITY 
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The general action for a f(R) theory of gravity is given 
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where R stands for the Ricci scalar, k := 8ttGq, and 
ip represents collectively the matter fields (we use units 
where c = 1). The field equation arising from varying 
the action Eq. (TTJ with respect to the metric is 

faRab - ^fg a b ~ (V Q V 6 - g ab U) f R = nT ab , (2) 

where f R :— d R f and □ = g ah V a Vf,. It is straightfor- 
ward to write the above equation in the following way 
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where (Vi?) 2 := g ab (V a R){V b R). Taking the trace of 
this equation yields 
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where T :— T a a . Finally, using Eq. (g]) in ©, we find 
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Equations (j4} and ((5|) are the basic equations for f{R) 
gravity that we propose to treat in every application, 
instead of transforming them to STT. Now, several im- 
portant remarks are in order regarding this set of equa- 
tions. First, aside from the case f{R) = R, where the 
field equations reduce to those of GR, one is to consider 
functions f(R) such that f R , f RR > (i. e. monoton- 
ically growing and convex f(R) functions) in order to 
avoid potential blowups in the field equations. These 
two conditions have been considered previously; the first 
one leads to a a positive definite gravitational "constant" 
G e ff = Go/ f R when regarding f(R) theories as effective 
theories. The condition fp R > was suggested to avoid 
gravitational instabilities [l2|. Second, Eq. ([5]) supplies 
a second order equation for the metric provided that the 
Ricci scalar is considered as an independent field. This 
is possible thanks to Eq. Q which provides the informa- 
tion needed to solve Eq. ([5]) for the metric alone (given 
a matter source). 

As will become evident below, this approach is rather 
clean and free of the pathologies that can arise when 
mapping f(R) gravity to STT. 



III. STATIC AND SPHERICALLY SYMMETRIC 
SPACETIMES 

In order to give some insight to the usefulness of 
our approach, we consider a static and spherically sym- 
metric (SSS) spacetime with a metric given by ds 2 = 
—n{r)dt 2 + m{r)dr 2 + r 2 (dO 2 + sin 2 Odip 2 ), where the 
metric coefficients n and m are functions of the coor- 
dinate r solely. The Ricci scalar as well as the matter 



variables will also be functions of r. The analysis of this 
kind of spacetimes is interesting in many ways; however, 
we will focus here only on the issue about the existence 
(or absence thereof) of well behaved compact objects. 
Equation ^ yields 
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(where ' := d/dr). From the t — t, r — r, and 6 - 
components of Eq. <j5j> and using also Eq. © we find 
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Note that Eqs. and © are not independent. In fact, 
one has the freedom of using one or the other (see Sec. IIVI 
below). Now, from the usual expression of R in terms of 
the Christoffel symbols one obtains 

4n 2 m(m — 1) + mm! (An + rn') 



R = 



2r 2 n 2 m 2 
—2rnm(2n + rn") + r 2 mn' 2 



(10) 



As one can check by a direct calculation, using 
Eqs. d?])-© in Eq. ((TDJ leads to an identity R = R. This 
result confirms the consistency of our equations. When 
defining the first order variables Q n = n' and Q R := R', 
Eqs. (f6j) — have the form dy l /dr = J- l (r, y z ), where 
y l = (m,n,Q n , R,Q R ), and therefore can be solved nu- 
merically. As far as we are aware, the system ((S]) — © 
has not been considered previously (see Ref. [l]| for an 
alternative approach). These equations can be used to 
tackle several aspects of SSS spacetimes in f(R) grav- 
ity. Finally, we observe that for f(R) = R the above 
equations also reduce to the well known equations of GR 
for SSS spacetimes. When dealing with extended ob- 
jects, notably those described by perfect fluids, the above 
equations are supplemented by the matter conservation 
equations. The Bianchi identities imply the conservation 
of the effective energy-momentum tensor [which corre- 
sponds to the right-hand-side (rhs) of Eq. which to- 
gether with Eqs. ([3]) and {5{ lead to the equation of con- 
servation V a T ab = for the matter alone. So for a perfect 
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fluid with T a b = {p+p)u a Ub+g a bP, the conservation equa- 
tion leads to p' = —{p+p)n'/2n [where n' is given explic- 
itly by the rhs of Eq. ©]. This is the modified Tolman- 
Oppenheimer-Volkoff equation of hydrostatic equilibrium 
to be complemented by an EOS. 

In order to solve the differential equations some bound- 
ary conditions must be supplied, which in this case are 
rather regularity and asymptotic conditions. Regular- 
ity (smoothness) at r — implies the following ex- 
pansion near r = 0: <f)(r) — (f>o + 0i^ 2 /2 + C(r 4 ) 
(where <fi stands for m, n, or R). This implies m' = 
= n' = R' at r = 0. We choose m(0) = 1 (local- 
flatness condition) and n(0) = 1. The coefficients cf>o 
and <fii associated with m, n, R, and the matter vari- 
ables (and which correspond to the values of these quan- 
tities and its second derivatives evaluated at r = 0, 
respectively) are related to each other. For instance, 
from the above power expansion and from Eqs. ([§]), © 
and mn, we find n"(0) = /°-2/S^-4fcr°+i8fc:r» and 

R (0) = — — iftn , where the quantities at the rhs 

J J RR ._. 

are evaluated at r = Q. 

Now, as concerns the asymptotic conditions, we are in- 
terested in finding solutions that are asymptotically de 
Sitter, since the solutions are supposed to match a cos- 
mological solution that represents the observed Universe. 
Therefore, we demand that asymptotically R — > Ri, 
where Ri is a critical point (maximum or minimum) of a 
potential which is defined below. The value Ri allows to 
define the effective cosmological constant as A cff = i?i/4 
(like in GR with A)0. We mention that at the end of 
the numerical integration, wc rescale n(r) in order to get 
the canonical asymptotic behavior n(r) ~ 1 — A c ffr 2 /3. 
This rescaling amounts to redefine the t coordinate. The 
asymptotic behavior R — > Ri depends on the value of 
i?(0). The details of the correlation between i?(0) and 
Ri depend on the matter model, and once this latter is 
fixed (for instance, an incompressible fluid with a given 
central pressure) the value R(0) can be found by a shoot- 
ing method [14| . Since we are interested in finding an 
exterior solution of Eq. (|6]) with R — > Ri (and R' — > 0) 
asymptotically, we observe that such a solution exists 
if Ri corresponds to a critical point of the "potential" 
V{R) = -Rf(R)/3 + J R f(x)dx. That is, R 1 is a point 
where dV(R)/dR = (2/ — fnR) /3 vanishes [assuming 
Irr(Ri) 7^ in Eq. ©] . This potential is radically 



1 Note that in GR the field equations imply the algebraic relation- 
ship R. = —kT, which determine the value of R in terms of the 
matter content solely. In f(R) theories this relationship is dif- 
ferential [cf. Eqs. J4) and {6)]. and therefore the conditions for 
R(0), R"(0) and even for n"(0) are not fixed in advance by the 
matter content only. We can recover the usual expressions of GR 
when taking f(R) = R. 

2 One can easily show that the system Eqs. — (J9j in vacuum 
has the exact de Sitter solution n(r) = m(r) — 1 = 1 — A c ff r 2 /3, 
R = R 1 = const, with A cff = i?i/4 and Ri = 2f(Ri)/f a (Ri), 
which corresponds precisely where dV{R)/dR = 0, as defined in 
the main text. 



different from the scalar-field potential that arises under 
the STT map. Furthermore, V(R) is as well defined as 
the function f(R) itself. 

A technical difficulty that one faces when integrating 
the equations for neutron star models that asymptoti- 
cally match realistic values of the cosmological constant 
is that a fortiori two completely different scales are in- 
volved. On one hand p(0) ~ p n uc ~ 10 14 gcm -3 ; on the 
other A - 10" 29 gcm" 3 (A = A/G ). That is, there are 
around 43 orders of magnitude between the typical den- 
sity within a neutron star and the average density of the 
Universe. This ratio between densities naturally appears 
in the equations since the parameters which define the 
specific f(R) theory are of the order of A, while the ap- 
propriate dimension within neutron stars is /O nuc - So in 
units of Pnuc the cosmological constant turns out to be 
ridiculously small, while in units of A, p(r) and p(r) turn 
out to be ridiculously large within the neutron star. So 
far, the authors that have studied relativistic objects in 
f(R) gravity [D-d, Q have faced this kind of technical 
problem, and, in order to avoid it, they have only con- 
structed objects which are far from representing neutron 
stars embedded in a realistic de Sitter background. This 
means that either the background is realistic while the 
objects are ridiculously large (several orders of magni- 
tude larger than a real neutron star) and not very dense, 
or the objects are realistic but A e ff is far too large. Such 
objects are relativistic in the sense that their pressure is 
of the same order of magnitude as their energy-density, 
and the ratio between a suitably defined mass and radius 
GqM/TL is not far from 4/9. 



IV. NUMERICAL RESULTS 

In order to test our method, we used first 
f(R) = R — aR* In (1 + R/R*) (a,R* are positive 
constants; R* sets the scale), which was proposed in 
Ref. [9(. In that analysis the authors mapped the theory 
to the STT counterpart. However, unlike the Starobinsky 
model (see below), in this case the resulting scalar- field 
potential turns to be single valued and the authors of 
9] did not find any singularity within the object. Under 
our approach, we associate to this f(R) the potential 

V(R) = 2j* | (1 + R) (R + 6a -1)- 2a(3 + 2R)\n(l + R)\, 

where R = R/R*. This potential has several critical 
points. Figure [1] (top panel, right inset) depicts the 
potential for a = 1.2. Equations ©— © were integrated 
using a four-order Runge-Kutta algorithm by assuming 
an incompressible fluid with the same density and 
central pressure as in Q. By implementing a shooting 
method [lij], we found a solution for R that goes from 
the global minimum V(R\) of the potential to a positive 
value. The minimum at R\ corresponds to the de 
Sitter value and gives rise to the effective cosmological 
constant A e g = i?i/4. Figure Q] (top panel) depicts 
the solution for i?, the pressure (top panel, left inset) 
and the metric potentials (bottom panel) plotted up to 
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the value where n(r) reaches the cosmological horizon 
given by n(r^) = 0. We checked that asymptotically 
the metric potentials matched perfectly well the de 
Sitter solution n(r) — m(r) -1 = 1 — A c ffr 2 /3. We also 
verified that our solution corresponds to the one found 
in Q and also confirmed that no singularity whatsoever 
was encountered. Furthermore, we checked that the 
numerical solutions found using the system Eqs. ©-© 
and the system Eqs. ©,© an d © gave the same results 
within a relative error no larger than ~ 10~ 6 and that for 
both systems the identity Eq. (fTU)) was satisfied within 
a relative error of ~ 10~ 10 . In this way, we ensured 
that no mistakes were introduced in the code during 
the typing of the equations. We emphasize that for this 
model the condition fjm > is satisfied by construction 
but the condition /r > is not in general. However, 
for the solutions we found, this latter is always satisfied, 
particularly at R\. Having analyzed the solutions for 
this f{R), we turned our attention to the Starobinsky 

model given by /(#) = #- A#* jl - [l + (i?/i?*) 21 " 

with (3 — 1. This is the model analyzed in 
using the transformation to STT. For such a model, 
the associated scalar-field potential is multiple-valued 
(cf. Fig. 1 in H| and Fig. 2 in [4|). Despite the 
potential being ill-defined, the authors in [J] argue 
that in the particular region where the scalar-field is 
evaluated, the potential is nevertheless well defined. It 
seems suspicious that in d, @ the odd features of the 
potential are not mentioned at all. Even if one tried 
to ignore the region where the potential is pathological 
and consider only the single-valued region, there is no 
guarantee that the ill-defined region will not play any 
role in more general settings than the very specific case 
of static and spherically symmetric compact objects. 
Modifications of the Starobinsky model by adding a 
quadratic term do not cure this pathology, as one 
can see in Fig. 5 of JTa. As we stressed before, the 
solutions found in [4j-l6| gave rise to a controversy 
since in Q the authors found a singularity in the Ricci 
scalar within the extended object, while in 0, Q the 
authors did not. Following our approach, we get V(R) = 
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R - 4A - 2A 1 + R 



3A arctan(i?) 



This potential has a rich structure depending on the 
value of A (such structure arises also for the potential 
discussed above for different values of a) . Figure [2] (top 
panel, left inset) depicts the potential for three values 
of A. There is a value A cr ;t = 8/\/27 below which the 
potential has only one critical point (the global mini- 
mum) located at R = 0. In this regime (A < A cr j t ), we 
have found only oscillatory solutions around the global 
minimum. For A > A cr ;t, three critical points appear, 
with a minimum always located at R — (Rq). In view 
of this structure, several kinds of solutions are possible 
for any given value of A > A cr it - We have found solutions 
where R asymptotically approaches a minimum of V(R) 
at Ri ^ 0, the local maximum, and also solutions where 
R oscillates asymptotically around Rq. An asymptotic 



oscillatory behavior around Rq seems a priori not ade- 
quate to produce the required de Sitter background; we 
plan to analyze those solutions in more detail elsewhere. 
Here, we are only interested in showing a solution for 
which R goes to R\ (a local minimum). Figure [5] shows 
the numerical solution which is asymptotically de Sitter. 
As one can appreciate, no singularity in the Ricci scalar 
was encountered in this solution and R never crosses 
the two real- valued zeros of fjm (which correspond to 
R = ±i?*/v3) nor those of Jr at which blowups in the 
equations can be produced. Even if for the static and 
spherically symmetric solutions, one can avoid the zeros 
°f Jr an d /rr, the Starobinsky model does not satisfy 
the conditions /ij, Jrr > in general, and so, in other 
settings (e. g. cosmology) one has to keep this in mind 
when solving the equations. 
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FIG. 1. (Top) Ricci scalar and pressure (left inset) as func- 
tions of r taking p = po — const. — 5 x 10 7 Ri/Go, p(0) = 
0.3po, and Ri/R, ~ 1.405 as in The stars depict the 
location of the object's surface. Potential V(R) (right inset) 
in units of Rl, for the model f(R) = R - ai?„ln (1 + R/rJ) 
with a = 1.2. (Bottom) metric potentials associated with the 
solution for R. 



V. CONCLUSIONS 

We have devised a straightforward and robust ap- 
proach to treat f(R) metric gravity without resorting to 
the usual mapping to scalar-tensor theories. With this 
method, it is possible to analyze in a rather transparent 
and well defined way several aspects of these alterna- 
tive theories. In particular, we focused on the existence 
of relativistic extended objects embedded in a de Sitter 
background and concluded that for some f(R) models 
such objects can be constructed without ambiguity and 
without resorting to any dubious explanations based on 
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FIG. 2. Similar as Fig. [T] for the Starobinsky model, with 
A = 1.56 and /3 = 1 taking p = 10 6 i?i/(167rG ), p(0) = O.lpo 
with .Ri/T?* ~ 1.983. 



the use of ill-defined quantities. It seems that the anal- 
ysis of the solutions presented here as well as other (i 
e. solutions that admit different de Sitter backgrounds, 
even possibly Minkowski backgrounds) can only be car- 
ried out easily following our method since our potential 
V(R) depicts many features in a rather clear and clean 
fashion that allows us to identify the critical points which 
R can reach asymptotically. 

Building realistic neutron stars with a realistic de Sit- 
ter background in f(R) gravity still remains a technical 
challenge. In the near future, we plan to study in more 
detail this issue along with other aspects of f(R) gravity 
using our approach. 
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